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Abstract 

We present a preliminary measurement of the time-dependent CP asymmetry for the neutral B- 
meson decay into the CP = +1 final state i^gTr'^vr'^, with vr+vr^ and vr*^ 77. We use 

a sample of approximately 227 million i?-meson pairs recorded at the T(4S') resonance by the 
BABAR detector at the PEP-II B-Factory at SLAC. From a maximum likelihood fit we extract the 
mixing-induced CP-violation parameter S'^o^o^o = 0.84 it 0.71 (stat) it 0.08 (syst) and the direct 
CP-violation parameter C^o^OttO = 0.27 ± 0.52 (stat) it 0.13 (syst), where the first uncertainty is 
statistical and the second systematic. 
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1 INTRODUCTION 



CP violation effects in decays of B mesons that are dominated by 6 ^ sqq transitions, where 
q = u,d,s, are potentially sensitive to contributions from physics beyond the Standard Model 
pp. The i3-factory experiments have explored time-dependent CP-violating (CPV) asymmetries 
in several such decays H], including 5° ^ ^K^ H Ej, B° K^K^K'^ ^° ^ v'K^ El, 
B° K+K-K^ ISHU, -BO ^ /o(980)A'° "B! and B^ K°7r0 ^. Within the Standard Model the 
asymmetry in these decays is expected to be consistent with the asymmetry in 6 — > gs decays, such 
as B^ J/ipKg, where the CPV asymmetry occurs due to a phase difference between mixing and 
decay amplitudes. These comparisons must take into account contributions of other amplitudes 
with different weak-interaction phases within the Standard Model. A major goal of the S-factory 
experiments is to reduce the experimental uncertainties of these measurements and to add more 
decay modes in order to improve the sensitivity to beyond-the-Standard-Model effects. 

In this letter we present a preliminary measurement of the CPV asymmetry in the decay B^ 
KgTT^TT^, using data collected with the BABAR detector at the PEP-II asymmetric-energy e~^e~ 
collider. In the Standard Model this decay is dominated by the b — > sqq amplitude, with q = u,d. 
A possible contribution from a tree-level b — > uus amplitude is doubly Cabibbo-suppressed with 
respect to the leading gluonic penguin diagram. 

The i^OvrOvrO final state is a CP-even eigenstate, regardless of any resonant substructure In 
the Standard Model we expect Sj^o^o^o — — sin 2/3 and Cj^o^o^o — 0. The angle /3 is defined as /3 = 
arg{—VcdV*i^/VtdV^l) and Vij are the elements of the CKM matrix [2]. A significant measurement 
of CP violation in this channel alone in comparison to other penguin modes constrains certain 
extensions of the Standard Model ^21 ■ 

2 THE BABAR DETECTOR AND DATASET 

The data in this analysis were collected with the BABAR detector |13j at the PEP-II asymmetric 
e~^e~ collider jl4j . A sample of 226.6 it 2.5 million BB pairs was recorded at the T(4S') resonance 
(center-of-mass energy y/s = 10.58 GeV). The BABAR detector is described in detail elsewhere |13) . 
Charged particles are detected and their momenta measured by the combination of a silicon ver- 
tex tracker (SVT), consisting of five layers of double-sided detectors, and a 40-layer central drift 
chamber, both operating in the 1.5 T magnetic field of a solenoid. Charged-particle identification 
(PID) is provided by the average energy loss in the tracking devices and by an internally reflecting 
ring-imaging Cherenkov detector (DIRC) covering the central region. Photons and electrons are de- 
tected by an electromagnetic calorimeter composed of 6580 CsI(Tl) crystals; the typical resolution 
for the TT^ signal in the 77 invariant mass spectrum is better than 7 MeV/c^. 

3 ANALYSIS METHOD 

In the decay B^ — > i^^vr'^Tr'^, which has no charged tracks originating from the B^ decay vertex, 
we rely on the technique recently developed to reconstruct the B^ vertex in B^ K^tt^ decays 
(described in detail below) 9 . From a candidate BB pair we reconstruct a B^ decaying into the 
CP eigenstate Kgir^-K^ (Bcp)- We also reconstruct the vertex of the other B meson (Btag) and 
identify its flavor. The difference At = tcp — itag of the proper decay times is obtained from the 
measured distance between the Bcp and Btag decay vertices and from the boost {(3'y = 0.56) of the 
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e^e system. The At distribution is given by: 

e-|At|/r 

r±{At) = -^—[l^Aw± (1) 
(1 - 2w) {Ssm{AmdAt) - C cos(AmdAt))]. 

The upper (lower) sign denotes a decay accompanied by a (B^) tag, r is the mean B^ hfetime, 
Anid is the mixing frequency, and the mistag parameters w and Aw are the average and difference, 
respectively, of the probabilities that a true B^ is incorrectly tagged as a B^ or vice versa. The 
tagging algorithm |E] has seven mutually exclusive tagging categories of differing purities (including 
one for untagged events that we retain only for yield determinations). The analyzing power, defined 
as efficiency times (1 — 2w)'^ summed over all categories, is (30.5 it 0.6)%, as determined from a 
large sample of i?-decays to fully reconstructed flavor eigenstates (-Bflav)- 

We search for B^ K decays in BB candidate events selected using charged-particle 

multiplicity and event topology [TB]. We reconstruct — > 7r"'"7r~ candidates from pairs of oppo- 
sitely charged tracks. The two-track combinations must form a vertex with a probability greater 
than 0.001 and a vr^vr^ invariant mass within 11.2 MeV/c^ of the nominal mass jl7| . We form 
vr" — > 77 candidates from pairs of photon candidates in the EMC, each of which is isolated from 
any charged tracks, carries a minimum energy of 30 MeV, and has the expected lateral shower 
shape. Candidates for B^ — > K^ti^tt^ are formed from K^ti^t^^ combinations and constrained to 
originate from the e^e~ interaction point using a geometric fit. We require that the consistency 
of the fit, which has one degree of freedom, be greater than 0.001. We extract the decay length 
Lj^o and the vr" 77 invariant mass from this fit and require 110 < m^^ < 160 MeV/c^ and L^o 
greater than 5 times its uncertainty. The cosine of the angle between the direction of the decay 
photon in the center-of-mass system of the mother vr" and the tt^ flight direction must be less than 
0.92. 

We extract the signal yield, S and C from an unbinned extended maximum likelihood fit where 
we parameterize the distributions of several kinematic and topological variables for signal and 
background events in terms of probability density functions (PDFs). 

For each B candidate we compute two kinematic variables, the energy difference AE = E*Q — ^^/s 

and the beam-energy-substituted mass ttt-es = \l {\s + po ■ pbY /Eq — p'^ , where s is the center- 
of-mass energy squared. The subscripts and B refer to the initial T{4S) and the Bcp candidate, 
respectively, and the asterisk denotes the center-of-mass frame. For signal events, AE is expected 
to peak at zero and thes at the known B mass. From a detailed simulation we expect a signal 
resolution of about 3.6 MeV/c^ in niEs and 45 MeV in AE. Both distributions exhibit a low-side 
tail due to the response of the EMC to photons. We remove a small dependence of the signal AE 
resolution on the location in the K^tt^t:^ Dalitz plot by using AE/a{AE) instead of AE, where 
a{AE) is the measured uncertainty in AE. We select candidates with ruEs > 5.20 GeV/cP and 
—5 < AE/a{AE) < 2. To suppress other B decays we also require —0.25 < AE < 0.1 GeV, which 
does not affect the signal AE/a{AE) distribution. 

The background B candidates come primarily from random combinations of and neutral 
pions produced in events of the type e+e~ qq, where q = u,d,s,c (continuum). Background 
from BB events may occur either in charmless decays B^ — > K^X, or from decays where the 
is from an intermediate charmed particle. The shapes of event variable distributions are obtained 
from signal and background Monte Carlo (MC) samples and high statistics data control samples. 
In rriEs, the charmless B background exhibits a broad enhancement near the i?-meson mass while 
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other B background distributions show no peaking. In AE/a{AE), B backgrounds in general show 
no clustering. 

In continuum events, particles appear mostly in two jets. This topology can be characterized 
with several variables computed in the T(4S') frame. One such quantity is the angle 9t between the 
thrust axis of the Bcp candidate and the thrust axis formed from the other charged and neutral 
particles in the event, where the thrust axis is defined as the axis that maximizes the sum of the 
magnitudes of the longitudinal momenta. This angle is small for continuum events and uniformly 
distributed for true BB events. With the requirement | cos 9t\ < 0.9 we suppress background by 
a factor of three while retaining 90% of the signal. We also use the angle 9b between the Bqp 
momentum and the beam axis, and the sum of the momenta pi of the other charged and neutral 
particles in the event weighted by the Legendre polynomials LQ{9i) and L2{9i) where 9i is the 
angle between the momentum of particle i and the thrust axis of the Bcp candidate. We combine 
these three variables in a neural net (NN) that is trained and evaluated [I^ on different subsets 
of simulated signal and continuum events and on data taken about 40 MeV below the nominal 
center-of-mass energy. The A^A^ has two hidden layers with 4 neurons each. The A^A^ output is 
divided into 10 consecutive intervals, chosen such that they are uniformly populated by the signal 
events; the PDF is modeled as a parametric step function ^H] whose parameters are the heights 
of each bin. Since the parent distribution for the A^A^ output is unknown any assumed functional 
form will suffer a systematic uncertainty due to the choice of the function. 

We suppress background from other B decays by excluding several invariant mass intervals: 
m(i^°7r°) > 4.8 GeV/c^ eliminates B° i^°7r°, 1.75 < m(if°7r°) < 1.99 GeV/c^ reduces B° 
D^TT^ to fewer than 10 expected candidates, m(7r'^7r*') < 0.6 GeV/c^ removes rjK'^ and rj'Kg, and 
3.2 < m{7r^Tr^) < 3.5 GeV/c^ removes XcoK^ and Xc2-f^s candidates. 

From MC simulation we expect more than one candidate in 13% of the signal candidate events. 
Because the number of multiple Kg candidates is negligible (less than 0.1%), we select the candidate 
whose two reconstructed vr'^ masses are closest to the expected value. The signal reconstruction 
efficiency is about 15%. 

For each B^ i^°7r°7r° candidate we examine the remaining tracks in the event to determine 
the decay vertex position and the flavor of i?tag- We parameterize the performance of the tagging 
algorithm in a data sample (-Bfiav) of fully reconstructed B^ — > D^*^~tt^ / /af decays. For the 
continuum background, the fraction of events tagged in category k, e^, is extracted from a fit to 
the data. The -Btag vertex is reconstructed inclusively from the remaining charged particles in the 
event (16; . 

To reconstruct the Bcp vertex from the single trajectory we exploit the knowledge of the 
average interaction point (IP), which is determined every 10 minutes from the spatial distribution 
of vertices from two-track events. The uncertainty on the IP position, which follows from the size 
of the interaction region, is about 150 //m horizontally and 4 //m vertically. We compute At and 
its uncertainty from a geometric fit [201 to the T(4S') B^B^ system that takes this IP constraint 
into account. We further improve the sensitivity to At by constraining the sum of the two B decay 
times {tcp + ttag) to be equal to 2 r^o with an uncertainty of -v/2 r^o, which effectively constrains 
the two vertices to be near the T{4S) line of ffight. This procedure provides an unbiased estimate of 
At. The extraction of At with the IP-constrained fit has been extensively tested on large samples 
of simulated B^ K decays with different values of S and C, and in data [2]. 

The per-event estimate of the uncertainty on At reflects the strong dependence of the At 
resolution on the flight direction and on the number of SVT layers traversed by the decay 
daughters. In about 70% of the events both pion tracks are reconstructed from at least 4 SVT hits. 
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leading to sufficient resolution for the time-dependent measurement. The average At resolution 
in these events is about 1.0 ps. For events that fail this criterion or for which a{At) > 2.5 ps or 
At > 20 ps, the At information is not used. However, since C can also be extracted from flavor 
tagging information alone, these events still contribute to the measurement of C. 

By exploiting regions in data that are dominated by background, and simulated events for 
the signal, we have verified that with our selection the observables are sufficiently independent 
that we can construct the likelihood from the product of one-dimensional PDFs, apart from the 
signal TTT-Es and AE / a{AE) which are correlated away from their mean signal positions and for 
which we use a two-dimensional PDF derived from a smoothed, simulated distribution. We obtain 
the PDF for the At of signal events from the convolution of Eq.(^ with a resolution function 
lZ{5t = At — Attrue, <7At)- The resolution function is parameterized as the sum of two Gaussians 
with a width proportional to the reconstructed ciAt, and a third Gaussian with a fixed width 
of 8 ps [T^. The first two Gaussian distributions have a non-zero mean, proportional to a^t, to 
account for the charm decays on the Stag side. We have verified in simulation that the parameters of 
TZ{6t, (7 At) for — > i^'^vr^vr*^ events are similar to those obtained from the Sflav sample, even though 
the distributions of a At differ considerably. We therefore extract these parameters from a fit to the 
-Bflav sample. We use the same resolution function for background from other charmless B decays. 
The At distributions for background from B decays into open charm final states and continuum 
consist of a prompt component and a non-prompt component, and the resolution function has the 
same functional form as used for signal events. The parameters for the At PDF of the open-charm 
background are determined from MC simulation, while for the continuum they are varied in the fit 
to data. 



4 MAXIMUM LIKELIHOOD FIT 



We subdivide the data into the tagging categories k, events with and without At information 
(set I and II), and those located in the inside or outside region of the Dalitz plot (inside or 
outside). The latter accounts for the higher contribution and different characteristics of continuum 



background near the Dalitz plot boundary. We define the quantity S = min{m 



2 

12' 



m 



13' 



m23), where 



rriij is the invariant mass of the B decay daughters i and j combined. It corresponds to the 
distance of an event in the Dalitz plot to the nearest Dalitz plot boundary in the limit of massless 
daughters. We split the data at 6 = 3.5 GeV^/c^. We maximize the logarithm of the extended 
likelihood C = e^^s+^s) . ^j^}^ JV5 and Nb (= J2b ^b) the total signal and background yields, 
respectively. The likelihood in each tagging category k (with tagging fraction e^) is given as: 
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The probabilities {Q^) and [Q^) for each measurement j are the products of PDFs 
for signal (5) and background {B) classes: P^j = PDF{mESj, '^E/a{AE)j) ■ PDF{NNj) ■ 
PDF{Atj, cj(At)j , tagfcj, kj), where for the background PDF{mESj, AE/a{AE)j) = PDF{mESj)- 
PDF{AE/a{AE)j). The probabilities Q do not depend on At and a{At) and are used to extract 
C from the yields. The fractions of events with At information for signal and background are 
denoted by fg and fg , respectively, and fractions of events in the outside Dalitz plot region by 
and For about 22% of our signal B candidates one or two of the tt^ decay photons associated 
with Be p originate from the i?tag • According to Monte Carlo simulation studies in these cross-feed 
events we expect to measure the same S and C as in the correctly reconstructed signal (true) since 
the contribution of the vr*^ to the At measurement is marginal. To account for differences in the 
PDF distributions for the signal probabilities P'^ {Q'^) we use: P = fcfPcf + (1 — fcf)Ptrue- The 
fraction of cross-feed events, /c/, is fixed to the value obtained from the simulation. Parameters of 
signal PDFs are the same for the different Dalitz plot regions. The PDFs for B backgrounds are 
identical for the Dalitz inside and outside regions. The tagging fractions for the signal and the B 
decay backgrounds are the same; continuum background has different . The good fractions fg 
and and the outside fractions and for continuum are varied in the fit, while these frac- 
tions for charm and charmless B backgrounds are determined from Monte Carlo simulations. The 
fit was tested with both a parameterized simulation of a large number of data-sized experiments 
and a full detector simulation. 

5 PHYSICS RESULTS 

The central values of 5 and C were hidden until the analysis was complete. From a data sample 
of 33,058 S° K^tt'^tt^ candidates, we find Ns = 117 ± 27 signal decays with S^^o^o^o = 0.84 ± 
0.71 (stat) ± 0.08 (syst) and Cj^o^o^o = 0.27 ± 0.52 (stat) ± 0.13 (syst) . The linear correlation 
coefficient between the two CP parameters is 2%. The yield of charmless B background is consistent 
with zero. Figure^shows the distributions of the event variables mES) AE/a{AE), and A^A'^, and 
Fig. [21 shows the At distributions for the B^- and the -B^-tagged subsets with the raw asymmetry 
[Nbo - N^o]/[Nbo + A^^o]. The A^^o (A^^o) is the number of B^ (B^) -tagged events. In all plots 
data are displayed together with the result from the fit after applying a requirement on the ratio 
of signal likelihood to signal-plus-background likelihood (computed without the variable plotted) 
to reduce the background. 

6 SYSTEMATIC STUDIES 

We consider systematic uncertainties listed in Table ^ These include the uncertainties in the 
parameterization of PDFs for signal and backgrounds which were evaluated by varying parameters 
within one standard deviation or using alternative shape functions. The largest uncertainty for C 
is caused by the A^A^ shape for continuum inside the Dalitz plot (ct(C) = 0.10) and for S from the 
2-D parameterization (cr(5) = 0.04). We consider uncertainties in the background fractions and CP 
asymmetry in the charmless B background, the parameterization of the At resolution function and 
the vertex finding method, knowledge of the event-by-event beam spot position, imprecision in the 
SVT alignment, and the possible interference between the suppressed b ^ u4 amplitude with the 
favored b ^ l^d amplitude for tag-side 5-decays [2j|. We fix r^o = 1.532 ps and Anid = 0.505 ps~^ 
and vary them by one standard deviation ^| . We correct for the small fit bias which is determined 
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NN shape Ls/(Ls+Lb) 




Ls/(Ls+Lb) 

Figure 1: Distribution of the event variables (a) rriES; (b) AE/a{AE), and (c) NN output in 10 
bins after reconstruction and a requirement on the ratio of signal likelihood to the signal-plus- 
background likelihood, calculated without the plotted variable. The solid line represents the fit 
result for the total event yield and the dotted line for the total background. Plot (d) shows the 
ratio of the signal likelihood to signal-plus-background likelihood with all variables included, data 
(dots) with the fit result superimposed. Plot (e) shows the same quantity as (d) close to one and 
with a linear scale. 
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Figure 2: Plots (a) and (b) show the At distributions of 5°- and -B^-tagged — > K^tt^tt^ 
candidates. The sohd hues refer to the fit for ah events; the dashed hnes correspond to the total 
background. Plot (c) shows the raw asymmetry (see text). A requirement is applied on the event 
likelihood to suppress background. 



from repeated fits to simulated events for signal and backgrounds mixed together with the expected 
yields, and the uncertainty of the method is accounted for as systematic error. 

We perform several consistency checks, including the measurement of the lifetime; we obtain 
r^o = 1.25 lb 0.47 ps. We embed different B background samples from Monte-Carlo simulation in 
the data sample and obtain consistent yields and CP parameters from the fit. We use the PDFs 
to generate signal and background samples and find that 47% of the simulated experiments had 
likelihood values greater than the one obtained in the fit to the data. 
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Table 1: Sources of systematic uncertainty on S and C. The total error is obtained by summing 
the individual errors in quadrature. 



Source 


a{S) 


a{C) 


PDF parameterization for signal and background 


0.05 


0.11 


Background fractions 


0.03 


0.02 


CP in charmless B background 


0.03 


0.01 


Vertex finding/Resolution function 


0.02 


0.05 


Beam spot position 


0.00 


0.00 


SVT alignment 


0.02 


0.01 


Tag side interference 


0.00 


0.01 


Amrf, TB 


0.02 


0.01 


Fit Bias 


0.04 


0.02 


Total systematic error 


0.08 


0.13 



7 SUMMARY 

We have presented a preliminary measurement of the CP violating asymmetries in Kgir^-K^ 
{K'^ — > TT~^iT~) decays reconstructed from a sample of approximately 227 million BB pairs. From 
an unbinned extended maximum likelihood fit we obtain 5^o^o^o = 0.84 it 0.71 (stat) it 0.08 (syst) 
and Cj^o^Oj^o = 0.27 it 0.52 (stat) it 0.13 (syst). The change in the log- likelihood when we fix the 
values of —£'^0,^07^0 to the average sin2/3 measured in b ^qs modes, sin2/3 = 0.725 ±0.037 [22j, and 
Cj^o^o^o to zero, and re-fit the data sample is 2.5. The signal yield is consistent with our findings 
in the B^ — > Kg7:~^7r~ decay [231 assuming the dominant charmless final states are /o(980)-ftr^, 
K*{892)7:^, i^Q (1430)7r'^, and non-resonant K^tt^t:^, and isospin symmetry. 
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